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Abstract. If one is willing to give up the cherished hypothesis of spatial isotropy, many 
interesting cosmological models can be developed beyond the simple anisotropically expand¬ 
ing scenarios. One interesting possibility is presented by shear-free models in which the 
anisotropy emerges at the level of the curvature of the homogeneous spatial sections, whereas 
the expansion is dictated by a single scale factor. We show that such models represent viable 
alternatives to describe the large-scale structure of the inflationary universe, leading to a kine¬ 
matically equivalent Sachs-Wolfe effect. Through the definition of a complete set of spatial 
eigenfunctions we compute the two-point correlation function of scalar perturbations in these 
models. In addition, we show how such scenarios would modify the spectrum of the CMB 
assuming that the observations take place in a small patch of a universe with anisotropic 
curvature. 


Contents 


1 Introduction 1 

2 Geometry and spatial eigenfunctions 3 

2.1 Anisotropy through spatial curvature 3 

2.2 Spatial eigenfunctions 4 

2.2.1 Flat FLRW 4 

2.2.2 Bianchi III 5 

2.2.3 Kantowski-Sachs 6 

3 Scalar perturbations 7 

3.1 General formalism 7 

3.2 Inflation and metric perturbations 9 

4 Sachs-Wolfe effect and correlation functions 10 

4.1 FLRW case 10 

4.2 Bianchi III case 11 

4.3 Kantowski-Sachs case 13 

5 Qualitative analysis 15 

6 Conclusion 16 

A Mathematical complements 17 

A.l Summation theorems 18 

A.2 Approximating Legendre polynomials by Bessel functions 18 

References 19 


1 Introduction 

The combination of observational data from modern galaxy surveys with those from the 
Cosmic Microwave Background (CMB) radiation strongly suggests that, on scales around 
100 Mpc and above, the spatial distribution of the cosmic web is homogeneous and isotropic. 
This highly symmetric description of the universe comprises what is known as the cosmolog¬ 
ical principle, and is mathematically encoded in the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) spacetime metric. 

Besides being backed up by observational data, the use of the FLRW metric has a 
technical advantage. Since the only free parameter in this metric is a function of time, 
background dynamical equations become ordinary differential equations, which can be solved 
by means of known analytical or semi-analytical techniques. More important, though, is 
the fact that, in this description of the spacetime, “pictures” of the three-dimensional spatial 
section at different times are copies of one another, differing only by a global conformal 
factor. If one further expands cosmological observables in terms of the eigenfunctions of 
the spatial Laplace-Beltrami operator, the aforementioned simplification implies that the 
eigenmodes of physical quantities will not couple through the evolution of the universe, as far 
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as the perturbative regime is maintained, and then one can trace their dynamics individually 
by means of linear differential equations. The easiness and versatility of this program has 
enormously boosted the field of cosmology in the past decades, culminating with the so-called 
standard cosmological model (the ACDM model), in which the qualitative features of the 
large-scale universe can be understood by means of only six free parameters. 

As comfortable as this description might seem, it is important to keep in mind that the 
cosmological principle constitutes a set of working hypotheses, and, in reality, a particularly 
critical set, given the intrinsic human limitations to perform observations at cosmological 
scales. In other words, one might be careful to not let arguments of simplicity prevent the 
discovery of new cosmological phenomena. 

The necessity of scrutinizing the symmetry hypothesis of the standard cosmological 
model also has the important role of testing the robustness of cosmological observations, 
which could otherwise be biased by a posteriori assumptions. As a matter of fact, while the 
observed approximate isotropy of the CMB is in agreement with the cosmological principle, 
this does not mean that different symmetry principles are ruled out by it and, in fact, different 
principles are known to be compatible with the CMB observations. Indeed, in the past few 
years several authors have considered the possibility that the spectrum of CMB could be fitted 
by large void models in which we would be located near its center [1-4], in clear contrast to the 
homogeneity hypothesis. The spectrum of CMB fluctuations has also been contrasted with 
spatially homogeneous but anisotropic cosmological models, either as a mean of constraining 
the impact of spatial anisotropy on CMB data [5-7], or as an attempt to explain known 
statistical anomalies of the CMB [8-11]. Of these two different routes, the latter is the 
one which we plan to follow in this work. That is, we assume spatial homogeneity of the 
universe and investigate the impact of spatial anisotropies on the large-angle spectrum of 
CMB temperature fluctuations. 

The simplest example of an anisotropic spacetime corresponds to a generalization of 
the FLRW metric in which each spatial direction expands according to its own scale factor. 
Such geometries are known as Bianchi I metrics. Given that the CMB signal points to a 
great degree of spatial isotropy of the early universe, Bianchi I metrics can be deployed in 
two different ways. For example, one could suppose that the onset of inflation has taken 
place in an anisotropic spacetime, and ask for its impact on primordial quantum fluctuations 
re-entering the horizon at the time of CMB formation [12]. On a different context, the 
possibility of a late-time anisotropization of the spacetime could leave the CMB spectrum 
essentially unchanged, while affecting both the growth rate of structures [13] as well as weak- 
lensing observables [14, 15]. As in the case of FLRW geometries, Bianchi I metrics also enjoy 
the simplicity of having unperturbed quantities that obey ordinary differential equations. 
However, the anisotropic expansion of the spacetime inevitably leads to a coupling of Fourier 
modes, which renders the use of this metric quite ingenious at the technical level. 

However, once we are willing to abandon the isotropy hypothesis of the cosmological 
principle, much richer models can be constructed besides the Bianchi I solution [16-20]. In 
fact, it has been demonstrated in references [18, 21-23] that the anisotropy of the universe 
does not need to emerge from its dynamical expansion; instead, it could result from the 
curvature of spatial sections, while having the expansion controlled by one scale factor, at 
the cost of having an imperfect fluid sourcing the anisotropic curvature. This idea has been 
successfully applied in the context of late-time anisotropies in references [24-26]. From the 
phenomenological point of view, this possibility offers two main advantages. First, since 
the background metric is controlled by a single scale factor, the redshift of electromagnetic 
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radiation is exactly isotropic, which means that the CMB remains exactly isotropic at the 
background level. Second, the fact that there is a single scale factor implies that pictures 
of the universe at different times are also copies of each other, exactly as in the FLRW 
description. Then, in the perturbative regime, eigenmodes of physical quantities will also 
evolve independently of one another, something which greatly facilitates the computation of 
cosmological observables. 

In this work we consider two spacetime metrics presenting isotropic expansion but 
anisotropic spatial curvature. These are the metrics from the Bianchi type III and Kantowski- 
Sachs families. Following the formalism of linear and gauge-invariant perturbations intro¬ 
duced in [27], we evaluate the impact of these metrics on the large-angle spectrum of CMB 
fluctuations, assuming that the particle horizon during CMB formation is much smaller than 
the curvature scale of these spaces. After introducing an explicit set of spatial eigenfunctions 
in Section 2, we show in Section 3 that, provided that the universe has gone an early period 
of exponential expansion, the dynamics of scalar perturbations will be described by one single 
gravitational potential $. From this result, an important corollary follows: given that large- 
angle CMB fluctuations result from the geodesics of photons in a perturbed spacetime, the 
Sachs-Wolfe effect in the considered geometries will have exactly the same functional form 
as in the FLRW case. This is demonstrated in Section 4, where, using the explicit eigen¬ 
functions found in Section 2, we compute the temperature two-point correlation function in 
these spaces, as well as their perturbative series around the isotropic two-point correlation 
function in the limit of large curvature radius. The effect of anisotropic curvature leads to 
off-diagonal multipolar correlations in the spectrum of the CMB which extends to a large 
multipolar range. This signature is investigated numerically in Section 5. We conclude with 
some remarks and perspectives of further developments in Section 6. 


2 Geometry and spatial eigenfunctions 

We start this section by introducing the main mathematical and technical tools which will be 
needed later on to compute the angular power spectrum of the CMB. 


2.1 Anisotropy through spatial curvature 

In this work we are interested in spatially homogeneous but anisotropic solutions to Einstein 
equations which nonetheless display an isotropic expansion, i.e., solutions whose dynamics 
is described by a single scale factor. Rather, the anisotropy of these solutions arises from 
the spatial curvature of the sections of constant comoving time, curvature which is not the 
same in all directions. This can be achieved by considering spatial sections which are the 
Cartesian product of curved subspaces. Since one-dimensional manifolds are trivially flat, the 
only possibility left for a three-dimensional manifold to present anisotropic spatial curvature 
is to be the product of a two-dimensional and a one-dimensional manifolds. In the case 
of maximally symmetric two-dimensional manifolds, and again apart from trivial flat cases, 
we only have to consider the possibility of a two-sphere or a bi-dimensional pseudo-sphere. 
Therefore, the spatial topology turns out to be x M or x M, respectively, accepting the 
topology of the real line for the third spatial dimension. The metric of the spaces that we are 
going to consider is thus conformally (ultra-)static and can be written as 


ds^ = a^{r]) —dr]^ -|- jabdx^'dx^ + dz^ 
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where {a, b} represent coordinates in the two-dimensional subspace and z is the coordinate 
of the real line. For future reference it will prove convenient to rewrite the line element in 
cylindrical coordinates 


ds^ = a^{r]) [—dr/^ -|- dp^ -|- S‘^{p)d(p‘^ + dz^] , 
where the function S{p) is given by 


S\p) 


< sinh^ p , 
_ sin^ p. 


(M3) 

{If X M) 
(§2 X M) 


( 2 . 2 ) 


(2.3) 


and where, to facilitate comparison with standard results, we have included the flat FLRW 
case (M^) as well. Furthermore, note that while 2 ; G (— 00 ,-|-cx)) and p G [0, 27r) in all three 
cases, the variable p belongs to the interval [0, -|-oo) in the and xM cases, while p G [0, tt) 
for X M. Moreover, note that we can always recover the flat space limit by taking p <C 1, 
which gives 

S\p)^p\ 

In the cosmology literature, metrics described by either the x M or the x M spatial slices 
are known as Bianchi type III (Bill) and Kantowski-Sachs (KS) solutions, respectively. 


2.2 Spatial eigenfunctions 

In order to compute correlation functions between physical observables (for instance, the 
gravitational potential), it is very convenient to deal with their spatial dependence by ex¬ 
panding them in modes that decouple in the dynamical held equations. These held equations 
are typically of a generalized Klein-Gordon type. Hence, for the case of scalar observables, 
we will hrst need to explicitly hnd and normalize (with respect to the inner product given by 
the volume element of the spatial sections) a complete set of scalar modes which are solutions 
of the eigenvalue problem 

V^Q = -q^Q^ qeR (2.4) 

where Q is a scalar function, and 






h = det{hij), 


(2.5) 


is the Laplace-Beltrami operator in a three-dimensional Riemannian space with metric hij. 
The eigenfunctions in the three cases considered here can be found by means of a simple 
separation of variables 

Q{p,P,z) = f{p)g{p)h{z). 

Below, we present these normalized eigenfunctions. 


2.2.1 Flat FLRW 

In the hat FLRW case the spatial metric is hij = diag(l, p^, 1). Using separation of variables, 
the normalized solutions to the eigenvalue problem (2.4) are easily found to be 

^imip ^ikz 

Zumki^) = W^/^Jm(wp)—, (2.6) 

p Itt 2tt 
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where the eigenvalues are 


m € Tj , A: G M, 


w G 


and obey the dispersion relation 


z 2 I 7 2 
q = io + k 


The normalization constant is such that 




Any scalar function Q(x) can be expanded in terms of this basis as 

fOO /* + co 


with inverse given by 


poo ^+CO 

e(x)= E / do; / dkQmioj, 

Q^{u;,k) = J \//id^xQ(x)Z;^^fc(x). 


(2.7) 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


2.2.2 Bianchi III 

In the Bill case the spatial metric in cylindrical coordinates is hij = diag(l, sinh^ p, 1). The 
normalized solutions in this case are 


^ikz 

Zuimk{^) — l_|_j^(cOshp)—=—= , 

^J2^l ^/2-K 


( 2 . 11 ) 


where 

(jj G M+, m G N, fc G M, 


and 


+ 


1 

4 ■ 


( 2 . 12 ) 


Note that, owing to the curvature of the spatial sections, there is a strictly positive lower 
bound to the frequency of a mode: 


q 


2 



(2.13) 


Modes with 0 < < 1/4 do not belong to the basis of square-integrable eigenfunctions of 

the Laplace-Beltrami operator, and are thus excluded from the frequency spectrum. Such 
frequencies lead to what is known as supercurvature modes [28]. 

The Legendre function (cosh p) - sometimes also called conical function - is real 

and has the following symmetries 


P 




T{iu) + m + 1/2) 


{ioj — m + 1/2) 
The (real) normalization constant Nmuj is given by 


P 


(^>1) 


(2.14) 


N„ 


(—l)™'a;tanh vrw 


r(ia; — m + 1/2) 
r{icj + m + 1/2) ’ 


(2.15) 
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and is such that [29] 


I \/hd^xZ^mfc(x)^*/^/fc/(x) = - w')(5(A: - k'). (2.16) 

Since these eigenfunctions form a basis for the Hilbert space of square-integrable functions 
with the volume element corresponding to hij, any such function Q(x) can be expanded as 


Q(x) 




dkQ 

m (w, k) , 


(2.17) 


with the inverse following directly from the orthogonality relation (2.16). The flat-space limit 
is reached when p <C 1 and the frequency is large, w ^ 1, such that cop is kept fixed. Using 
the approximations 


(ia;)-/’-™/^(coshp) ~ i-irJrr^icop), 

Tjiuj + m + l/2) ^ 2 m 
r(it<; — m-|-1/2) ’ 

one can check that, up to an overall and unimportant phase, 

^im(p ^ikz 

^Lomki^) ^ UJ ^ Jm{uJp)—7==—7== , 

V ZTT V ItT 


(2.18) 

(2.19) 


( 2 . 20 ) 


as it should be. 


2.2.3 Kantowski-Sachs 

Finally, for the KS case, the spatial metric in cylindrical coordinates is hij = (l,sin^p, 1). 
Since the spatial sections are just the product of the unit two-sphere with the real line, the 
eigenfunctions are simply the standard spherical harmonics multiplied by a one-dimensional 
plane wave: 


Zemki^) = NmePPicosp) 


^im(f ^ikz 


where 


and we now have 


£gN, meZ, k eR, 
= £{£ + !) +k'^. 


( 2 . 21 ) 


( 2 . 22 ) 


Although the £ = 0 mode is included in the basis of square-integrable functions, this mode 
corresponds to an overall monopole, and is excluded from the definition of cosmological per¬ 
turbations. Thus, the lower bound to the frequency of a mode is 


q^>2. 


The normalization constant is defined as 


^ I i2£ + l)i£-m)\ 
y 47r(.^-|-m)! 


(2.23) 


(2.24) 
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and is chosen to give 


j \/hd^xZ£mfc(x)Z|,^,fc,(x) = 6w6mm'S{k - k'). (2.25) 

Since these solntions also form a complete basis for the space of sqnare-integrable fnnctions, 
we can again expand any snch fnnction as 

^ ^ r+oo 

Q(x) = ^ ^ / dkQtm{k)Zimk{?^) ■, (2.26) 

£=0 

with the inverse following from (2.25). The flat-space limit may again be achieved by letting 
p <C 1 and 1, with ip held constant. Using the identity = {—l)'^{i — rn)\/{i+m)\P'P‘ 
and the approximations P^~™'(cosp) r\j i ^Jm{ip) and (i + m)\f{i — m)\ ~ we find 

Akz 

Zimk{-^) ~ i^/^J^{ip)^^ , (2.27) 

V 27r V 27r 


again, np to a phase. 

3 Scalar perturbations 

In this section we recall the basic formalism of linear pertnrbation theory in Bill and KS 
backgronnds [27]. Then, by analyzing the linearized Einstein eqnations of the pertnrbative 
degrees of freedom we show that, provided that the nniverse has gone an early period of 
exponential expansion, scalar pertnrbations can be parameterized by a single free fnnction, 
exactly as in the standard FLRW case. Onr main resnlt is (3.16). Readers not interested 
in the details of the derivation can skip them and go directly to Section 4. Thronghont this 
section we nse the coordinates introdnced in formnla (2.1). 

3.1 General formalism 

The theory of linear cosmological pertnrbations is bnilt npon the idea that the large-scale 
nniverse can be described by means of small pertnrbations added to a backgronnd (exact) 
solntion of Einstein eqnations. In general terms, this program can be implemented with the 
introdnction of ten arbitrary fnnctions on top of the backgronnd spacetime metric. Then, 
for homogeneons and isotropic spatial sections, these fnnctions can be classified as either 
Scalars, Vectors or Tensors (SVT) by exploring their transformation properties nnder the 
action of the spatial symmetry gronp. In an FLRW nniverse this leads to the known SVT 
mode decomposition [30]. The main advantage of this procednre is that, in the linear regime, 
each kind of pertnrbation evolves independently of the others, which simplifies the analysis of 
their dynamical evolntion. In this decomposition, non-trivial SVT modes exist in all FLRW 
cosmologies whose spatial slices are copies of a maximally symmetric space. In the Bill or 
KS spacetimes, on the other hand, spatial slices are given by the prodnct of the real line with 
a two-dimensional maximally symmetric space. Therefore, in this case it seems natnral to 
classify the pertnrbations modes jnst according to their behavior with respect to the symmetry 
gronp of that two-dimensional manifold [27]. In this sense, it is worth remarking that there is 
no non-trivial transverse and traceless symmetric tensor in two dimensions. More specifically. 
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on the two-dimensional maximally symmetric subspaces, any vector I 4 and symmetric tensor 
Tab can be decomposed as 

Va = DaV + Va, (3.1) 

Tab = X'-Jab + DaDbY + D^^a^b) j (3-2) 

where Da is the covariant derivative compatible with the metric 'jab [see (2.1)] and barred 
vectors are transverse, i.e., DaV^ = 0 = DaZ°^. Let us also notice that, similarly, no non¬ 
trivial transverse vector can exist in the remaining one-dimensional section of the spatial slices, 
in which only scalar modes can be defined. Actually, this decomposition of the perturbations is 
an adaptation of the more general formalism of Gerlach and Sengupta [31], extended beyond 
spherical symmetry, where we treat the two dimensions supplementary to those of the bi- 
dimensional maximally symmetric manifold not in a covariant way as a Lorentz manifold on 
its own, but rather described in terms of the conformal time and the z-coordinate introduced 
above. 

Since the background metric of the Bill and KS models considered here have a single 
scale factor, there will be no dynamical coupling between SV modes, something which further 
simplifies the dynamical evolution of the system. As stated in the Introduction, this situation 
is in sharp contrast with perturbation theory in Bianchi I universes where, although one 
can define uncoupled SVT modes at a given initial time, the anisotropic evolution of the 
background metric couples those modes, leading to a dynamical see-saw effect [12, 32]. Once 
clarified that perturbative mode couplings will not arise in the present case, we can focus 
exclusively on the evolution of scalar perturbations, which are the most important ones for 
the CMB large-angle perturbations. The most general metric with scalar perturbations can 
be written 

ds^ = o^(r/)[—(1 -|- 2<h)dr7^ — 2Yi<lr]dz -|- (1 -|- 2T)7afedx“dx^ + (1 + 2A)d2;^] , (3.3) 

where <I>, 11, T, and A represent the four effective scalar degrees of freedom of the metric.^ 

Next, we need to introduce perturbations to the matter sector. Such perturbations can 
be parameterized in the same way as above, with the addition of extra SV degrees of freedom 
to the matter components. One important aspect of the formalism is that, since the considered 
models do not have a dynamical FLRW limit, the anisotropic spatial curvature needs to be 
sourced by an imperfect matter component [21]. This can be done either by introducing a 
stress component in the energy-momentum tensor [21] or a two-form field [24]. We can also 
take the simpler route provided by the choice of scalar fields, which can be chosen to be a 
real massless field in the Bill case [22, 23], or a complex massless field in the KS case [33]. 
Thus, our matter sector is composed of two fields: the inflaton, whose perturbations can be 
introduced exactly as in the FLRW case (see, for example, [34]), and the above anisotropic 
and massless scalar field <p(rj,x.). Their perturbations are as follows: 



(3.4) 


(3.5) 

a^6T^^^, = {A-dMK, 

(3.6) 


(3.7) 

~ + dzS4’ ) 

(3.8) 


^ Seven scalar functions in the original parameterization minus three scalar gauge functions [27]. 



where the prime denotes the derivative with respect to rj. As we are going to show, the 
pertnrbation 6(p of the anisotropic scalar held is not dynamically important in the inhationary 
stage, and can be safely ignored. 

In conclnsion, given the above pertnrbations to both the geometrical and the matter sec¬ 
tors, linear and gange-invariant dynamical eqnations can be straightforwardly compnted [27]. 

3.2 Inflation and metric perturbations 

We are primarily interested in compnting large-angle effects (> 1°) of the Bill and KS 
geometries on the temperatnre angnlar spectrnm of the CMB. At these scales, temperatnre 
pertnrbations are mainly driven by qnantnm hnctnations of the early nniverse - a period 
which, we assnme, was sonrced by scalar helds. Then, from the a — b component of the 
linearized Einstein eqnation, and from the fact that scalar helds do not prodnce anisotropic 
stress, it follows that (see (A.5) in [27]) 

- (A-h <I>) = 0 , a^b. (3.9) 

The above eqnation translates into the following constraint: 

A = -^>. (3.10) 

This is in close analogy to the sitnation of standard (FLRW) pertnrbation theory, where the 
absence of anisotropic stress is nsed to eliminate one scalar degree of freedom of the metric. 

Given the constraint (3.10), metric scalar pertnrbations are still described by three free 
fnnctions: <I>, T and H. Snrprisingly, thongh, only one degree of freedom is relevant dnring 
inhation. In order to see that, let ns consider the z — a component of Einstein eqnations: 

-{d,^ + d,^) + nu + ^u'= 6(1), (3.11) 

where Ti is the Hnbble parameter. If inhation is sonrced by a slowly-rolling scalar held, the 
variable 64> obeys [27] 


w" -F 


' 2 2 + 3e 

q + —^ 


w = 0. 


w = aScj ), 


(3.12) 


where e = l — is the slow-roll parameter and ry G (— oo, 0“) is the conformal time dnring 

inhation. Note that, since the “mass” term (2 -|- 3e)/ry^ is strictly positive, w will never grow, 
which in tnrn implies that Scf) will be snppressed by the expansion. Indeed, for ^ 77 “^, 
w ~ exp[±iqr]] and 6(1) decays as 1/a. For <C 7?“^, w ~ exp[±iy/3€ + 7/41og(—7y)]/(—ry)^/^, 
which implies that^ 6(1) ~ 1/a^/^. Since this resnlt holds for both Bill and KS, we can set 


6(f> = 0 


(3.13) 


with no loss of generality. 

If we now rewrite (3.11) in momentnm space (9^ —)• ik) and nse (3.13), we hnd that the 
constraint 

ikMq, rj) + ^{q, ry)] = (3.14) 

^We recall that, during slow-roll inflation, a ~ 1/{—(]) at dominant order in e. 
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must hold for all values of k. Moreover, from the dispersion relations (2.12) and (2.22) we see 
that q and k can be treated as independent variables whose specihc values hx uj (Bill) or i 
(KS). In other words, the case /c = 0 in the constraint (3.14) must hold for all possible values of 
n(g, r]). This is only possible if iq) is just a function of g, that is, if n(( 7 , rf) = /(g) /a?‘{ri). 

We thus see that the variable 11 is dynamically suppressed as well, and can be ignored. Since 
this result holds for both Bill and KS geometries, we can set 

n = 0, T = -4>. (3.15) 

We have thus arrived at the important conclusion that, provided that the universe 
has experienced an early period of exponential expansion, the metric describing large-scale 
perturbations in both Bill and KS geometries depends on a single free function, and can be 
parameterized as 

ds^ = a^(r/)[—(1 -|- 2<I>)dr/^ + (1 “ 2<I>)hjjdxMx'^] . (3.16) 

This is formally the same as the metric of adiabatic scalar perturbations in FLRW universes. 
However, our result is more general because hij{x^) can represent the FLRW, the Bill, or 
the KS spatial metrics. 

4 Sachs-Wolfe effect and correlation functions 

The main conclusion of the previous section is encoded in (3.16), and it implies that scalar 
perturbations over the considered manifolds with anisotropic spatial curvature are described 
by a single gravitational potential, exactly like standard perturbation theory in FLRW uni¬ 
verses. Moreover, (3.16) implies that all dynamical equations are formally the same as in the 
standard isotropic case - the only difference being that covariant derivatives are now taken 
with respect to the more general metric hij. A corollary of this result is that the Sachs-Wolfe 
(SW) effect, being purely kinematical, is given by the same known expression 

Ar(h) = ^4 >(x,77) , (4.1) 

where h is the direction of detection of a photon emitted at spacetime position (x, r/) (i.e., 
the coordinates of the last-scattering surface) and, again, x are the coordinates in the more 
general spacetimes under discussion. In this description, the effect of the anisotropic spatial 
curvature appears, in particular, through the non-trivial spectrum of the Laplace-Beltrami 
operator, and affects the angular power spectrum of the CMB at large scales, as we are now 
going to see. 

4.1 FLRW case 

We move now to the computation of the temperature two-point correlation function (2pcf) 
resulting from the Sachs-Wolfe effect in a flat FLRW universe. While this result is well-known, 
this section is intended to introduce the main formalism, and also to facilitate comparison 
with new results which will follow below. Moreover, here we will carry out the computation 
using cylindrical coordinates, that are not the ones most commonly employed. 

The temperature correlation function between two CMB photons arriving from the last 
scattering surface, considering only the SW effect, is dehned by 

(AT(h)Ar(fi')> = = ^(cl>(x)4.(x')), (4.2) 
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where, for simplicity, we have dropped the instants of emission, since these should be clear 
from the context. We now use (2.6) to decompose the gravitational potential ^>(x) and use 
the fact that, in a Gaussian and statistically isotropic universe, the correlation between modes 
takes the form 

k')) = - w')5(fc - k'). (4.3) 

Here 'P{q) is the primordial power spectrum and q is dehned as in (2.7). Employing the 
summation theorem for Bessel functions (A.l), we then arrive at 

1 poo p-\-oo 

C'(n, h') =-Tv-^ / wdw / dA:P(g) Jo(a;Apo)e**'^^ (4.4) 

(6vr)^ Jo 

where we have introduced the dehnitions 


(Apo)^ = + p‘^ ~ ‘^PP’ cos Ap , A(p = ip — Lp', IS.Z = z — /J . (4-5) 


As we will see, expressions (4.4) and (4.5) play a central role in the derivation of results 
for geometries with anisotropic curvature. Just for completeness, though, let us note that the 
right hand side of (4.4) can be further simplihed by the change of variables 


oj = q sin and k = q cos ^|J , 


(4.6) 


and the subsequent integration over the range^ 0 < tp < ir. The fact that the 2pcf depends 
only on the angle between h and h' can be seen by introducing spherical coordinates through 


= ArisinO^'^ and z^'^ = ArjcosO^'^ , 

where Ary is the conformal distance to the last scattering surface. This leads to 

1 


C(h,h') = C(r9) = 


ISvr^ 


2 J N [gAry\/2 - 2 cos r?] 
q dqV[q )'-- - 


'0 


qAr]\/2 — 2 cos r? 


(4.7) 


(4.8) 


where cos r? = cos 9 cos 6' + sin 9 sin 9' cos A<p. From here, the angular power spectrum Ci can 
be directly computed by means of a Legendre transform over r?: 


Ci = 27r / d(cosr?) Pi{cos'd) C{d) 


(4.9) 


1-1 


which, after performing the integral in cos r? by means of the change of variable 2y^ = 1 — cos r?, 
leads to the well known expression [30] 

= g^^(9) j|(9Ary) ding. (4.10) 

4.2 Bianchi III case 

In the Bill geometry the temperature 2pcf is still dehned by (4.2). However, the gravitational 
potential is now decomposed using the eigenfunctions (2.11). The perturbations in this case 
are supposed to be Gaussian but statistically anisotropic, so that the most general correlation 
function between modes will have the form 


(4'm(w, k)^*^,{u)\ k')) = P(q)hmm'<^(w - J)5{k - k '). (4.11) 

^The range of ip is conditioned by the restriction cj > 0. 
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Since we are working with homogeneous spaces, both Kronecker and Dirac deltas are needed 
to ensure translational invariance [35]. The power spectrum, on the other hand, is a general 
function of the vector q = {u},m,k), since the underlying geometries are anisotropic. In 
principle, a rigorous quantization procedure is required to fix and normalize 'P(q) properly. 
Nonetheless, note that invariance under rotations in the maximally symmetric bi-dimensional 
subspace requires a power spectrum independent of m. Moreover, since we are interested 
in departures from the CMB features typical of an FLRW universe resulting from the Bill 
geometry, it suffices to fix 'P(q) by demanding that the two-point correlation functions in 
these geometries agree in the limit in which the two points coincide (limit in which the 
spatial curvature is ignorable). This condition leads to 

tanhTTo; ’ 

where V is the same function appearing in (4.3) and where we have defined = 

( 7 ^ — 1/4 in order to account for the discrepancy between the dispersion relation in the FLRW 
case and (2.12). The factor tanhvra; is needed to ensure the correct functional dependence 
of the 2pcf when h —)• n/ according to our arguments. Using the above expression and the 
summation theorem (A.3) for Legendre polynomials, we arrive at 


iP(q) 




1 

(Gvr)^ 




d/cP(g)P_i_^.^(cosh A/3)e*^^^ , 


(4.13) 


where 

cosh Ap = cosh p cosh p' — sinh p sinh p cos A(p . (4-14) 

Before continuing our analysis, it is interesting to compare these formulas with their 
counterparts in the FLRW geometry. First, since P_i/ 2 +i(.j(l) = 1 = To(0), we can see 
that (4.13) agrees with (4.4) when Ap = Apo = 0 = Az, as they should. The main differences 
between these expressions are: (i) the law of cosines, equation (4.5), which now becomes the 
corresponding law in a two-dimensional hyperbolic space, namely (4.14), and (ii) the Bessel 
function in the kernel of the integral (4.4), which is now replaced by a Legendre function. 
As we will see, all expressions agree at first order in the limit where Ap <C 1 in units of 
curvature scale (and hence, strictly speaking, even beyond the limit of coincident points that 
we commented above). 

Apart from our physical requirements of coincident points, expression (4.13) is still quite 
general, inasmuch as it provides the form of the correlation between any two scalar functions 
in a space with x M topology. However, we are only interested in situations in which this 
2pcf does not lead to large modifications in the CMB data from the results of the flat FLRW 
case, since we know from observations that our universe is well fitted by those data. With 
this motivation, we can think of the CMB observations as taking place in a region described 
just by a small patch of a larger Bill region. If that is the case, small curvature corrections 
from the global Bill metric might be lurking beyond the particle horizon scale, and their 
signatures should be encoded in small corrections to the large-angle CMB spectrum. Such 
corrections come from two different terms. The first of them is the function T’_ i^j^(cosh Ap), 
which can be written as a powers series in Ap (see Appendix A. 2): 


-P_i+ia;(coshAp) = Jo{uAp) 


(M! 

12 


Jo{uAp) 


JijujAp) 

2u)Ap 


(4.15) 
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The second correction conies from the law of cosines (4.14), which can be expanded in powers 
of p and p' as 

(Ap)2 = (Ap„)2 + 55^(w')" + '" (4-16) 

where Apo was dehned in equation (4.5). Both expansions above are exact, and can be carried 
at any desired order in their perturbation parameters. In our case, this parameter is the 
distance to the last scattering surface in units of the curvature scale, Arj, which relates to Ap, 
p and p' through equations (4.7), (4.5) and (4.16). Notice that this is the only free parameter 
appearing in our equations. If we thus combine equations (4.15) and (4.16), retaining only 
corrections of order (Ap)^, we arrive at 


P_i+i^{coshAp) K. Jo(x) + 


[pp'Y sim A(p 


dJo{x) 

Apo? 

Jo(x) 

Ji{x) 

dx 

12 

2x 


6(Apo)2 

where x = a;Apo- Inserting this result in equation (4.13) we hnally have that 


, (4.17) 


C(n,fi') 


r+oo 


ikAz 


(6vr) 


+ 


(67r)^ 


wdcj / dkV{q)Jo{ijjApo)e 

J — OO 


6(Apo)^ Jo 


d{ujApo) 




Jo(wApo) — 


JijojApo) 

2ujApo 


AkAz 


. (4.18) 


Expression (4.18) is one of our main results. Note that the effect of the anisotropic 
curvature will appear both at the isotropic (h = h') and anisotropic (h ^ h') parts of the 
CMB spectrum. We will investigate these signatures with more details in Section 5. 


4.3 Kantowski-Sachs case 

The computation of the 2pcf in the KS case follows the same lines as for Bill. Starting 
from (4.2), we decompose the gravitational potential using the eigenfunctions (2.21). Since 
we are considering a homogeneous but spatially anisotropic geometry, the correlation between 
modes must take the form 


= VAAiAmm'Ak - k') , (4.19) 

where q = (£, m, k). Using arguments similar to those given in the previous section, one can 
show that 

:P(q) = 27rP(g). (4.20) 

Here we have dehned if = k"^ + {I + 1/2)^ = (7^ + 1/4 to account for the difference between 
the FLRW dispersion relation and (2.22), and we have chosen the numerical factor in the 
power spectrum so as to formally reproduce the expression for the 2pcf of the hat FLRW 
case in the limit of coincident points. Using the above power spectrum and the summation 
theorem (A.5) we hnd 


C'(n,n ) = 


(dvr) 


E 

i=i 




r+oo 


dkV{q)Pi{cos Ap)e 


ikAz 


(4.21) 
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where 


cos A/? = cos p cos p' + sin p sin p' cos{p> — ip). 


(4.22) 


A remark is in order, concerning the derivation of (4.20) and (4.21). In the limit of coincident 
points, equation (4.4) becomes (67r)~^ f coduj f dkV(q). Clearly, the integral over uj in this 
expression can be split as a sum of all the integrals over two consecutive positive integers. 
Assuming that the function V does not vary much with lu in any such interval, we can 
evaluate it at its middle point and then, in the considered case of coincident points, we are 
simply left with the integrals of the remaining factor wdcu. In the interval + 1], e.g., 
this integral is equal to + 1/2). Hence, with our assumption, the isotropic 2pcf becomes 
(bvr)”^ ^^(^ + 1/2) f dkP(q). The requirement that the 2pcf in the KS space should give 
an equivalent result in the limit of coincident points then justify (4.20). As a hnal comment, 
note that the main difference of the 2pcf in the KS space, apart from the modihcations in 
the dispersion relation and the eigenfunctions, is the appearance of a sum over £ instead of 
an integral over uj. This is reminiscent of the fact that the KS geometries have closed spatial 
subspaces in their spatial sections, leading to discrete eigenvalues of the Laplace-Beltrami 
operator. 

Moving forward, we follow the same line of reasoning as in the previous section and 
consider the relevant region for the CMB observations as a small patch of a larger KS space. 
Consequently, we expand P£{cos Ap) as a power series in Ap [see (A.9) in the Appendix] to 
hnd 


Pg{cos Ap) 


Jq{£Ap) + 


(Ap)^ 

12 


JoiiAp) 


JijlAp) 

2£Ap 


(4.23) 


where we have used the notation £ = £ + 1/2. Note that, in contrast to (4.15), the hrst 
non-trivial correction now appears with a positive sign. Likewise, we expand (Ap)^ dehned 
in (4.22) as a power series, which gives 


(Ap)^ 


(Apo)" - 


sin^ Ap 


{ppf + • • • 


(4.24) 


Combining these two expressions and retaining powers up to (Ap)^, as before, we hnd 


Pi{cosAp) Ri Jo(x) 


(pp')^ sin^ Ap 

6(Apo)2 


dJo{x) 

dx 


, (Apo)^ 

12 


Jo{x) 


Ji{x) 

2x 


(4.25) 


where x = {£ + l/2)ApQ. It is worth noticing that, since we are assuming that Apo oc Arj is 
small, the Bessel functions in these expressions do not vary much in each interval [£,£ + 1]. 
Thanks to this fact, our discussion below (4.22) still applies and allows us to regard the sum 
over £ as a discretized approximation of the integral over dcu. Combining all these results, we 
hnally hnd 


C(n,fi') 


r+oo 


(6vr) 


ojduj / d/c 7^(q)Jo(a;Apo)e 


ikAz 


(pp')^ sim Ap 


(67r)2 [ 6(Apo)2 

(Apo)2 


ujdcj dk V{q) (wApo) 


/o 


12 


/o 


J —oo 

oo ^+oo 

ujduj I dkV{q) 


d{ujApo) 


Jo(wApo) — 


JijojApo) 

2wApo 


AkAz 


, (4.26) 


with q^ evaluated at + A in these integrals. This is another of our main results. 
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5 Qualitative analysis 


We will now turn to a numerical analysis in order to extract some qualitative information from 
formulas (4.18) and (4.26). We start by noticing that these expressions can be summarized 
in a single formula as 

67(11, n') = (7(1?) ± J^(h, h'), (5.1) 

where C(i?), with •d = arccos(n • h'), is the isotropic 2pcf and where the anisotropic func¬ 
tion J^(h, h') is defined by the additional terms in the last two lines of equation (4.18) or 
equation (4.26). Here, the plus and minus signs correspond to the Bill and KS cases, re¬ 
spectively. Alternatively, our results can be written in terms of the temperature covariance 
matrix^ defined in terms of the temperature harmonic coefficients, a^m = f d^hAT(h)Y^'^(h), 
as 

= y J d^n C(h, h')y£m(h)y£/m'(n') • (5.2) 

Decomposing each term in C{h, h') into spherical harmonics, we find 

{Ofim^l'rn') ~ C£5££i6jnm' ^ ■ (5-3) 


There are two types of corrections coming from the anisotropy of the spatial curvature. 
The first of them arises from the diagonal part of the matrix {J-£m£'m’) ■ Such correction leads 
to an effective temperature power spectrum which differs from C£ mainly at small multipoles. 
Unfortunately this effect does not offer a strong constraint on the model since, besides being 
strongly limited by cosmic variance, it is degenerate with isotropic physical mechanisms which 
also alter the low-6 tail of the spectrum (e.g. the Integrated Sachs-Wolfe effect). The second 
type of correction,on the other hand, arises from off-diagonal terms in the covariance matrix, 
and is a genuine signature of statistical anisotropy. Moreover, this correction is present at 
all multipolar scales, a feature which alleviates the problem of cosmic variance. The simplest 
off-diagonal term we could expect to find is of the form ~ ■ Since 

this term couples even and odd multipoles, it is a signature of parity breaking. However, 
these terms are exactly zero in our case, since the metrics we are considering are invariant 
under the parity transformation (p,(p,z) —)• (/?, (/? ± tt,—z). The next off-diagonal term we 
can calculate is (otma(£+2)m)’ "''^^ere So, for each 6 there will 26-1-1 terms of 

this type. Recall also that a£m = i — as required by the reality of Ar(h). Then, 

to further reduce the effect of cosmic variance, we will work with an azimuthally averaged 
measure of anisotropy defined as a sum over the 26-1-1 possible realizations of the label m. 
Therefore, we define the following (real) quantity: 


J^e+2 = 


1 

26-6 1 


£ 

\{^irna{£+2)m)\ ' 

m=—i 


(5.4) 


We have evaluated this quantity numerically using a power spectrum per logarithmic band 
of the Harrison-Zel’dovich form: q^'P{q) = A = constant, and two different choices of the 
parameter Ar] compatible with our hypothesis that Arj <C 1. We show in figure (1) a plot of 
6(6-6 \)J -£-\.2 for a wide range of the multipole 6. Our analysis shows that the spectrum grows 
smoothly with 6 up to 6 = 100, where the use of the Sachs-Wolfe effect alone is justified. 
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Figure 1. Log-log plot of the quantity as a function of the angular multipole and in 

units of the power spectrum amplitude A. The curves show the results of our numerical computations 
assuming A?7 eqnal to 0.04 (dashed, magenta) and 0.05 (continuous, blue), in units of the cnrvatnre 
scale. 

6 Conclusion 

The hypothesis that we live in a spatially inhnite, homogeneous, and isotropic universe is 
one of the central pillars of modern cosmology. Besides being strongly supported by CMB 
measurements, the simplihcations introduced by this hypothesis in the framework of general 
relativity has allowed cosmologists to build a robust description of the large-scale universe with 
only six free parameters. Nonetheless, the fact that computations are simpler in the framework 
of an FLRW metric should not prevent us from considering more general geometries, especially 
when one realizes that the observed isotropy of the CMB does not rule out less symmetric 
models of the universe. Indeed, it is known that the anisotropies of the background metric can 
manifest themselves through other channels different from the simple anisotropic expansion 
of the scale factor. 

In this work we have explored the ideas introduced in references [21, 23, 24], in which the 
anisotropy of the universe does not result from its dynamical expansion, but rather from the 
curvature of its spatial sections. In the context of three-dimensional spaces, this feature can be 
implemented naturally only in two different ways: either by considering the Bianchi type III 
or the Kantowski-Sachs metric. In both models the expansion of the universe is controlled by 
a single scale factor, and hence the resulting CMB held is exactly isotropic at the background 
level [23]. Moreover, the fact that the dynamics are dictated by a single scale factor introduces 
simplihcations which render the model tractable from a computational point of view. Based 
on this property, and starting from the full theory of linear and gauge-invariant perturbations 
introduced in [27], we have established several important results. First, by means of a simple 
parameterization of the considered metrics in cylindrical coordinates, we have given explicit 
expressions for a complete set of spatial eigenfunctions possessing the symmetries of these 
spaces. This step is crucial for the comparison of the theory with observations and, to the best 
of our knowledge, has not appeared before in the context of anisotropic cosmologies. Second, 
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we have shown that, provided that the early universe has expanded in accordance with the 
inflationary paradigm, large-scale and adiabatic scalar perturbations will be described by one 
single gravitational potential, $, just as in standard perturbation theory in isotropic universes. 
Since the main large-angle CMB effects depend only on the photon geodesics in the perturbed 
spacetime, our result shows that the temperature fluctuations at these scales will have exactly 
the same functional dependence on $ as the one found in the isotropic setup. In particular, 
the Sachs-Wolfe effect in both Bill and KS spaces is given by AT = <h/3. 

Focusing on the Sachs-Wolfe effect, which is the most important effect for CMB at large 
angles, and with the help of the basis eigenfunctions that we have found, we have explicitly 
computed the spatial two-point correlation function in the Bill and KS geometries in the 
limit of small Ar], that is, a small ratio of the distance to the last scattering surface to the 
curvature radius. There are two different corrections to the isotropic correlation function, one 
resulting from the expansion of the spatial eigenfunctions around their isotropic counterpart, 
and another resulting from the law of cosines in curved spaces. The general effect of these 
corrections is two-fold: first, it changes the low-£ tail of the isotropic spectrum, and second, 
it induces off-diagonal temperature correlations which ought not to be present if the universe 
is isotropic. We have evaluated the latter numerically, and we have found that the effect goes 
beyond the low multipole region of the spectrum, corresponding to small i, region where, at 
least in principle, these correlations could be detected. It is important to note that we cannot 
predict the ratio of the last scattering surface to the curvature radius (i.e., Ar/), which in this 
work was treated as a small free parameter. Nonetheless, current bounds on super horizon 
features are usually of the order of a few percent [36] or even less [37], which means that our 
choices for Arj are based on educated guesses. 

Finally, it is worth noticing that the formalism that we have developed here is quite 
general and can be directly applied to situations which include other effects, such as the 
integrated Sachs-Wolfe and Doppler effects, as well as in scenarios involving more general 
matter fields. This work also paves the way to more challenging investigations, such as the 
effect of anisotropic curvature in the spectrum of gravitational waves and the polarization of 
the CMB. Such analyses, together with a thorough assessment of the detectability of these 
effects with actual and upcoming experimental probes, are currently in progress. 

Note Added After this work was submitted to publication we became aware of refer¬ 
ences [38-40] , in which the spatial eigenfunctions of B3 and KS spacetimes were found inde¬ 
pendently. We thank Julian Adamek for bringing these references to our attention. 
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A Mathematical complements 

In this appendix we succinctly review the main mathematical tools used in this work. 
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A.l Summation theorems 

The derivation of the 2pcf in flat (FLRW) geometries makes nse of the following identity 

+O0 

^ = Jo{oj/^Po), (A.l) 

m=—oo 


where Apo is defined in (4.5). The derivation of the 2pcf in the open (Bill) case, on the other 
hand, relies on the identities 


p—m 


^ ^ T{v + m + l) ^ 


+ 0O 

P_.+,„(co5hAp) = ^ _^(coshp)F™,^.„{cosh/)e‘’"<^-«’'> 

m=—oo 


(A.2) 

(A.3) 


where Ap is defined in (4.14). Analogonsly, the derivation of the 2pcf in the closed case nses 
the following formnlas: 

p,-” w = (-ir|(^pr w. (A.4) 

H-oo 

P^(cosAp)= ^ (-l)'"p-™(cosp)P^(cosp')e™^''■‘^'^ (A.5) 


where Ap is defined now in (4.22). 


A.2 Approximating Legendre polynomials by Bessel functions 

In order to find an approximation to P-i/ 2 +iujicosh.y) in terms of Bessel fnnctions we start 
with the following integral representations [41]: 


Jn{^) 


cos Lot 


V2 n 

vr Jo \/cosh y — cosh t 

{x/2r 


dt 


r(n + l/2)r(l/2) J_, 


1 - 


, y > 0, 

^/^cosxtdt. 


(A.6) 

(A.7) 


We now define t = uy, x = ooy, and write the argnment of integral (A.6) as a power series in 
y. This gives 


^/2 


"i 


y cos XU 


P_i+i^{coshy) ^ Ccoshy - coshny 


1 


r-l 


= — / COS XU 


TT 


'0 


du 

2/2(1 + ^2) 


= M^) - 


Lv^ 

Jo(x) - 


I2C1 - n 2 


+ 


du 


Jijx) 

2x 


+ • 


(A.8) 


Note that this series expansion is exact. In particnlar, since the resnlting series is controlled 
by y, we can trnncate it at any desired order in this variable if we regard a: as a fixed nnmber. 
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Following the same procedure as above, but now using 


Pi{cos9) - 
one can show that, defining x 


cos(£+l/ 2 )t 
TT Jo VCOS t — COS 9 
■ + 1/2)0, one has [41] 


0 < 0 < vr, 


P£{cos 0 ) 


q2 

Mx) + TW 


Jo{x) 


Jljx) 

2x 


(A.9) 


Again, note that the series is controlled by 0, for constant x. 

The reader familiar with the so-called flat sky approximation (usually employed in CMB 
analysis) will recognize expansion (A.9). It says that, for small patches of the CMB sky, the 
eigenfunctions of the sphere, i.e., the Legendre functions, converge towards eigenfunctions 
of the plane, i.e, Bessel functions [42]. For the same reason, expression (A. 8 ) gives the 
eigenfunctions of the pseudo-sphere in terms of eigenfunctions of the plane. 
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